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Abstract 

Vector-boson fusion processes constitute an important class of reactions at hadron 
colliders, both for signals and backgrounds of new physics in the electroweak interac- 
tions. We consider what is commonly referred to as W + W~ production via vector- 
boson fusion (with subsequent leptonic decay of the Ws), or, more precisely, e + v e yTv^ 
+ 2 jets production in proton-proton scattering, with all resonant and non-resonant 
Feynman diagrams and spin correlations of the final-state leptons included, in the 
phase-space regions which are dominated by i-channel electroweak-boson exchange. 
We compute the next-to- leading order QCD corrections to this process, at order a 6 a s . 
The QCD corrections are modest, changing total cross sections by less than 10%. 
Remaining scale uncertainties are below 2%. A fully- flexible next-to- leading order 
partonic Monte Carlo program allows to demonstrate these features for cross sections 
within typical vector-boson-fusion acceptance cuts. Modest corrections are also found 
for distributions. 



1 Introduction 



Vector-boson fusion (VBF) processes form a particularly interesting class of scattering 
events from which one hopes to gain insight into the dynamics of electroweak symmetry 
breaking. The most prominent example is Higgs boson production, that is the process 
qq — > qqH, which can be viewed as quark scattering via t-channel exchange of a weak boson, 
with the Higgs boson radiated off the W or Z propagator. Alternatively, one may view this 
process as two weak bosons fusing to form the Higgs boson. Higgs boson production via VBF 
has been studied intensively as a tool for Higgs boson discovery [TJ |2] and the measurement 
of Higgs boson couplings [3] in pp collisions at the CERN Large Hadron Collider (LHC). The 
two scattered quarks in a VBF process are usually visible as forward jets and greatly help 
to distinguish these Hjj events from backgrounds. 

An important background to Higgs searches at the LHC, in particular to the search 
for H — > W + W~ decays in VBF production, is caused by continuum W + W~ production in 
VBF. The qq — > qq W + W~ process forms an irreducible background in Higgs searches which 
ranges between 15% and 3.5% of the Higgs signal, for Higgs boson masses between 115 and 
160 GeV jl]. In fact, the kinematic distributions of the two tagging jets, the suppression of 
gluon radiation in the central region (due to the t-channel color-singlet exchange nature of the 
VBF process) and many features of the leptonic final state are identical to the H — > W + W~ 
signal. When trying to determine Higgs boson couplings, the qq — *> qqW + W~ cross section 
must be known precisely, which is achieved by calculating the next-to-leading order (NLO) 
QCD corrections. Such a calculation becomes more crucial when one contemplates using 
weak-boson scattering processes, and, more precisely, the absence of strong enhancements 
in these cross sections, as a probe for the existence of a light Higgs boson jJJ E|. Here the 
knowledge of NLO QCD corrections is essential in order to distinguish the enhancement from 
strong weak-boson scattering from possible enhancements due to higher order QCD effects. 

In two recent papers, the calculation of the NLO QCD corrections was presented for 
two simpler vector-boson-fusion processes: the Hjj signal cross section [7] and the cross 
sections for Zjj and Wjj production jSj. Both calculations were turned into fully- flexible 
parton-level Monte Carlo programs. We here extend this work and describe the calculation 
and first results for the NLO QCD corrections to W + W~jj production via VBF. 

Weak-boson scattering was first considered in the framework of the effective W approxi- 
mation, where the incoming weak bosons are treated as on-shell particles [D]. This approxi- 
mation does not provide a reliable prediction for the kinematical distributions of the forward 
and backward jets which are the main characteristic of vector boson fusion processes [TU] . 
Calculations of the full qq — > qq W + W~ processes, first without W decay [TTJ [TJ] and then 
including the full spin correlations of the W decay products in the narrow-width approxima- 
tion ^3], have been available for more than a decade. Within this latter approximation, also 
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the real gluon emission contributions, i.e. the 0(a A a s ) cross sections for the pp — > W + W~jjj 
subprocesses, with full spin correlations of the W decay leptons, were determined Jl] . Very 
recently, a partonic-level Monte Carlo for all the processes q±q 2 — > qzq±%%lv, with exact 
matrix elements at 0(a 6 ), has become available |To] . 

In this paper, we consider the proton-proton scattering process pp — > e + u e [i~v^ X , 
with all resonant and non-resonant Feynman diagrams and spin correlations of the final-state 
leptons included, at order a 6 a s . Since this process is very difficult to detect above QCD 
backgrounds, except in phase-space regions which are completely dominated by t-channel 
electroweak (EW) boson exchange, we only consider t-channel contributions, as explained 
in Sec. 12.11 In the rest of the paper, we will refer also to this approximated process as EW 
W + W~jj production. Electroweak gauge invariance requires that, beyond vector-boson 
scattering graphs, also the direct emission of the produced (virtual) Ws off the quark lines 
be considered. Several examples are depicted in Fig. [TJ which shows the basic Feynman-graph 
topologies which need to be considered for our calculation at tree level, for the particular 
subprocess uc~^uce + v e yTVp. Real emission contributions (including quark-gluon initiated 
subprocesses) are generated by attaching an external gluon in all possible ways on the two 
quark lines in Fig. ^ For the virtual corrections, we only need to consider Feynman graphs 
with a virtual gluon attached to a single quark line: gluon exchange between the up- and 
the charm-quark line leads to a color-octet state for the external uu or cc pair, which cannot 
interfere with the color-singlet structure at tree level. As a result, the virtual contributions 
contain, at most, pentagon diagrams, which arise e.g. by connecting the incoming and the 
outgoing up-quark in Fig. [T] (a) with a virtual gluon. The other graphs in Fig. [T] lead to 
box, vertex, or quark self-energy corrections, and these latter classes have already been 
encountered in Ref. 

Many aspects of the present calculation parallel this previous work. The cancellation of 
collinear and soft divergences for generic VBF processes was described in detail in Ref. [7j 
and need not be repeated here, since it can be applied verbatim for the case at hand. The 
calculation of vertex and box corrections was needed for the case of W and Z production [S] 
already, and, thus, these aspects of the virtual corrections need a brief review only. This 
review is provided in Sec. where we describe the details of our calculation and the approx- 
imations with regard to crossed diagrams in the presence of identical quark flavors. As in 
the previous work, we regularize the loop integrals via dimensional reduction and separate 
the virtual amplitudes into 1/e 2 and 1/e terms, which multiply the Born amplitude, and 
remaining finite terms, which are then calculated numerically, using the helicity-amplitude 
techniques of Ref. [T^j. A major concern here is the numerically stable and fast evalua- 
tion of the pentagon graphs. We make use of Ward identities and map large fractions of 
the pentagon contributions onto more easily evaluable four-point functions. Another im- 
portant feature is the systematic use of "leptonic tensors" which describe groups of purely 
electroweak subdiagrams. 
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In Section El we describe the numerous consistency tests which we have performed, 
ranging from comparison to code generated by MadGraph J7| for the tree-level amplitudes 
to gauge invariance tests. In addition, we present the properties of our numerical Monte 
Carlo program and how we have dealt with the gauge invariant handling of finite W/Z 
widths, the singularities for incoming photons and the choice of physical parameters. We 
then use this Monte Carlo program to produce first results for EW W + W~jj production at 
the LHC. Of particular concern is the scale dependence of the NLO results, which provides 
an estimate for the residual theoretical error of our cross-section calculations. We discuss 
the scale dependence and the size of the radiative corrections for various distributions in 
Sec. |U Conclusions are given in Sec. 03 

2 Elements of the calculation 

Our goal is the calculation of EW W + W~jj production cross sections with NLO QCD 
accuracy in phase-space regions which are typical for vector-boson fusion. This implies 
that some electroweak contributions, like triple gauge boson production (pp—>W + W~V with 
V—*jj), can safely be neglected. These approximations will be specified below. Also, we 
make use of the general structure of NLO QCD corrections to VBF processes: it is sufficient 
to specify the contributions to the Born, the real-radiation and virtual amplitudes which 
enter the cross section expressions of Ref. [7j. In this section, we describe how these contri- 
butions have been computed, the approximations used throughout this calculation and some 
technical details. 

2.1 Tree-level contribution and approximations 

The Feynman diagrams contributing to pp — > jj e + v e /z - ^, where both resonant and non- 
resonant processes are fully considered, can be grouped into six separate classes. The first 
group of two, which consists of the VBF processes considered in this paper, is characterized 
by t-channel neutral-current (NC) and charged-current (CC) exchange between the two 
scattering quark lines. The other four classes correspond to u- and s-channel exchange. The 
NC and CC labels are assigned depending on the external quark flavors: the incoming and 
outgoing quark charges on each quark line coincide for a neutral current process and differ 
by one unit of \e\ for a charged current process. 

For each neutral-current process, and in the unitary gauge which we use throughout, 
there are 181 Feynman graphs, which can be grouped into six distinct topologies. Generic 
diagrams for each of the topologies (a) to (f) are shown in Fig. [I] for the specific subprocess 
uc^uce + u e fi~v^. They correspond to the following configurations: 
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(e) (f) 

Figure 1: The six Feynman- graph topologies contributing to the Born process 
uc^uce + u e fi~u^, a template for neutral- current processes. Diagrams analogous 
to (a), (b), (e) and (f), with vector-boson emission off the lower quark line, are not 
shown. 



(a) Two virtual W bosons are emitted from the same quark line and in turn decay lepton- 
ically. 

(b) A virtual 7 or Z boson (V) with subsequent leptonic decay is emitted from either 
quark line. The tree-level expression for the sub-amplitude V — > e + u e is given by 
the tensor Ty, where a is the tensor index carried by the vector boson. 
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(c) The leptonically-decaying W bosons are emitted from two different quark lines. 



(d) Vector-boson fusion in the t-channel gives rise to the sub-amplitude VV — > e + u e /U - ^, 
which is characterized by the tensor Tyy. The tensor indices of the scattering V bosons 
are indicated with a and f3. 



(e) The leptons are produced by an external W boson emitted from a quark line and a 
W + V — ► e + i> e fusion process in the t-channel. The latter is described by T^f +V . 

(f) The leptons stem from W + emission from a quark line, accompanied by t-channel 
W~V — ► ii~v ll scattering, described by T^-v 



The propagator factors l/{q 2 — m v + imyTy) are included in the definitions of the sub- 
amplitudes introduced above, which we call "leptonic tensors" in the following. 




Figure 2: Diagrams contributing to the scattering amplitude T^ +v , depicted in 
Fig-U\ (e), that describe the tree-level subprocess W + V— >e + v e , where V is a 7 or 
a Z vector boson, and a and (3 are the tensor indices carried by the charged and 
neutral vector bosons, respectively. 



The explicit structure of one of these leptonic tensors is given in Fig. |21 where we have 
plotted the Feynman diagrams contributing to T^ +v : a virtual W + and a virtual 7 or Z 
fuse into a final state e + u e lepton pair, and the sub- amplitude corresponding to these three 
graphs is the leptonic tensor Tff+v wn i cn appears in graphs like Fig. [T] (e). 

For each charged-current process, such as us — > dce + i> e \i~Vp. or dc^us e + v e pTv^ there 
are 92 Feynman graphs. The different topologies are completely analogous to the ones for 
neutral current processes: simply interchange the t-channel bosons 7, Z *-> W in Fig. 
The only new tensor structure that occurs is T^ +w -, which describes the sub-amplitude for 
W + W~ — > e + u e n^u^. The corresponding Feynman graph topology is depicted in Fig. El 

By crossing the external quark lines, one either obtains anti-quark initiated t-channel 
processes like uc— >wce + z/ e yTv^ (which we fully take into account in our calculation) or one 
arrives at NC or CC s- or w-channel exchange between the two quark lines, which we count 
as the other four classes of ije + v e yU~i> M processes: 
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Figure 3: Contribution from W + W~ fusion to the scattering process 
us — > dce + u e pTVp. The tensor T^ +w _ contains all the tree-level diagrams con- 
tributing to the process W + W~ — > e + v e where a and (3 are the tensor indices 
carried by the W + and W~ vector bosons, respectively. 



- s-channel exchange leads to diagrams where all the virtual vector bosons are time-like. 
They correspond to diagrams called conversion, Abelian and non-Abelian annihilation 
in Ref. JH], and contain vector-boson production with subsequent decay into pairs of 
fermions. 

- u-channel exchange occurs for diagrams obtained by interchange of identical initial- or 
final-state (anti)quarks, such as in the uu—^uue + h , e fi~i , fl subprocess. 

In our calculation, we have neglected contributions from s-channel exchange completely. In 
addition, any interference effects of t-channel and w-channel diagrams are neglected. This is 
justified because, in the phase-space region where VBF can be observed experimentally, with 
widely-separated quark jets of very large invariant mass, the neglected terms are strongly 
suppressed by large momentum transfer in one or more weak-boson propagators. Color 
suppression further reduces any interference terms. In Ref. fE] we have checked that, for 
the analogous process pp^W/Zjj, the contribution from the two neglected classes and 
from interference effects accounts for less than 0.3% of the total cross section, at leading 
order. Since we expect QCD corrections to the neglected terms to be modest, the above 
approximations are fully justified within the accuracy of our NLO calculation. 



2.2 Real corrections 

The real-emission corrections to EW W + W~jj production with a gluon in the final state 
are obtained by attaching one gluon to the quark lines in all possible ways. There are 836 
graphs in the case of neutral-current processes and 444 for the charged-current ones. 

The contributions with an initial-state gluon are obtained by crossing the previous dia- 
grams, promoting the final-state gluon as incoming parton, and an initial-state (anti-)quark 
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as final-state particle. We again remove all diagrams where all electroweak boson propagators 
are time-like. Such diagrams, for consistency, must be removed since we have not consid- 
ered the corresponding Born contributions, namely the s-channel diagrams corresponding to 
triple weak-boson production. These diagrams are strongly suppressed when VBF cuts (see 
Sec. HJ) are applied to the final-state jets. 

In the regions of phase space where soft and collinear configurations can occur, we en- 
counter singularities in the phase-space integrals of the real-emission squared amplitudes. 
The regularization of these singularities in the dimensional-regularization scheme, with space- 
time dimension d = 4 — 2e, and the counter-terms which are needed to get finite expressions 
within the subtraction method, are discussed extensively in the literature (see, for exam- 
ple, [19J). Since these divergences only depend on the color structure of the external par- 
tons, the subtraction terms encountered for EW W + W~jj production are identical in form 
to those found for Higgs boson production in VBF [Zj and for EW Vjj production [Hj. The 
integration over the singular counter-terms yields, after factorization of the parton distribu- 
tion function, the contribution 



Here, the notation of Ref. [T2|, but adapted to dimensional reduction, has been used. M.b 
denotes the amplitude of the corresponding Born process and Q 2 is the momentum transfer 
between the initial and final state quark in Fig. ^ These singular terms are eventually 
cancelled by the virtual corrections, when infrared-safe quantities are computed. 

2.3 Virtual corrections 

As for the real-radiation cross sections, the divergences that affect the virtual gluon contri- 
butions depend on the color structure of the external partons. The main difference with Hjj 
and Vjj production is that the finite parts of the virtual corrections are more complicated 
for the present case, since the previous two processes only sport vertex and box corrections, 
while now we have to deal with pentagon-type loop integrals. 

The QCD corrections to EW W + W~jj production appear as two gauge- invariant subsets, 
corresponding to gluon emission and reabsorption on either the upper or the lower fermion 
line in Fig. ^ Due to the color-singlet nature of the exchanged electroweak bosons, any 
interference terms of the Born amplitude with virtual sub-amplitudes with gluons attached 
to both the upper and the lower quark lines vanish identically at order a s . Hence, it is 
sufficient to consider radiative corrections to a single quark line only, which we here take as 
the upper one. Corrections to the lower fermion line are an exact copy. We have regularized 
the virtual corrections in the dimensional reduction scheme (201 : we have performed the 
Passarino-Veltman (PV) |2Ij reduction of the tensor integrals in d = 4 — 2e dimensions, while 



< 1(e) >= \M B 
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the algebra of the Dirac gamma matrices, of the external momenta and of the polarization 
vectors has been performed in d — 4 dimensions. 



We split the virtual corrections into three classes: the virtual corrections along a quark 
line with only one vector boson attached (e.g. diagram (d) in Fig. [T]or diagrams (a), (b), (e) 
and (f) when considering corrections to the lower quark line), the virtual corrections along a 
quark line with two vector bosons attached (e.g. diagrams (b), (c), (e), (f)), and the virtual 
corrections along a quark line with three vector bosons attached (e.g. diagram (a)). 

I. The virtual NLO QCD contribution to any tree- level Feynman amplitude which 
has a single electroweak boson V\ (of momentum q x ) attached to a quark line, 

q{h)^q{k 2 )+V x {q x ), (2.2) 

is factorizable in terms of the amplitude for the corresponding Born graph 



M y = M u^ c (j^\ r(1 



2 3 

o - - + Cvirt + O (e) 



e 2 e 



(2.3) 



Here fiR is the renormalization scale, and the boson virtuality Q 2 = — {k\ — k 2 ) 2 is the only 
relevant scale in the process, since the quarks are assumed to be massless, k\ = k\ = 0. In 
dimensional reduction, the finite contribution c v i rt is equal to n 2 /3 — 7 (c v i rt = 7r 2 /3 — 8 in 
conventional dimensional regularization) . 

II. The virtual QCD corrections to the Feynman graphs, where two electroweak bosons V\ 
and V 2 (of outgoing momenta q x and q 2 ) are attached to a quark line, are depicted in Fig. HJ 
It suffices to consider one of the two possible permutations of V\ and V 2 , with kinematics 

q(h) -> q(k 2 ) + V 1 {q 1 ) + V 2 (q 2 ) . (2.4) 

Due to the trivial color structure of the tree-level diagram, the divergent part (soft and 
collinear singularities) of the sum of the four diagrams in Fig. 0] is a multiple of the corre- 
sponding Feynman graph at Born level, just like for the vertex corrections, 



47r V Q 1 



e z e 



+ ^-C F M^ lV2)T (Qi,q2) e 2 g^g? h + 0(e). (2.5) 

where we define Q 2 = 2k\ ■ k 2 , in order to use the same notation as in Eq. (|2.3|) . Here r 
denotes the quark chirality and the electroweak couplings eft* follow the notation of Ref . [TH] , 
with, e.g., g]f = Qf, the fermion electric charge in units of \e\, g^^ = 1 /{y/2 sin 9w) and 
gZf _ ^ jgy. — Qf sin 2 9 W ) / (sin Q w cos 9 W ) , where 9w is the weak mixing angle and T 3 f is the 
third component of the isospin of the (left-handed) fermions. 
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Figure 4: Virtual corrections for a fermion line with two attached vector bosons, 
Vi(qi) and V 2 {q2)- The finite part of the sum of these graphs defines the reduced 
amplitude M.YiVi,t iflli Q.2) of Eq. \2. 5)) . 



A finite contribution of the virtual diagrams, which is proportional to the Born amplitude 
(the c v i rt term) , is pulled out in correspondence with Eq. (j!2.3j) . The remaining non- universal 
term, M«y T (q±, q 2 ), is also finite and can be expressed in terms of the finite parts of the 
Passarino-Veltman .By, Cij and functions. The corresponding analytic expressions were 
given in Ref. [S]. Note that the effective polarization vectors for the electroweak bosons 
V\ and V~2, which enter the expressions for the M.y\v 2 r (Q11Q2), are W /± decay currents, the 
leptonic tensors Ty (for Fig. C](b)) and/or the entire lower parts of the Feynman graphs for 
Fig. E (c,e,f), when combining Feynman graphs with identical topology. 

III. The virtual QCD corrections to the Feynman graphs where three electroweak bosons 
Vi, V2 and V3 (of outgoing momenta qi, q 2 and ^3) are attached to a quark line, are depicted 
in Fig. El It suffices to consider one of the six possible permutations of V\, V 2 and V3, with 
kinematics 

q(h) -> q(k 2 ) + V 1 (q 1 ) + V 2 (q 2 ) + V 3 (q 3 ) . (2.6) 

The trivial color structure of the tree-level diagram allows the factorization of the divergent 
part of the sum of the eight diagrams in Fig. |3] in terms of the corresponding Born sub- 
amplitude 



o "T Cyirt 

e z e 



10 



k 2 k! 



qi| <i2£ q3? qi| q2| qs? qi| qi| q3. 
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qi? q2? q3? qi? qi? qs, 

Vj V 2 V 3 Vj v 2 v 3 

(g) GO 

Figure 5: Virtual corrections for a fermion line with three attached vector bosons, 
Vi(li)> ^2(92) & n d Vs^qs). The finite part of the sum of these graphs defines the 
reduced amplitude M-ViV2V 3 ,t (91,92,93) of Eq. j2.7\ ). 



asifiR) C F M {t) 



An 



V x V%Vz,t 



h, 92, 9s) e 3 gX xh 9? h 9X* h + 



(2.7) 



Again, a finite contribution from the virtual diagrams, proportional to the Born amplitude 



(cyirt), is pulled out and the remaining finite part is indicated with M. 





ViViVs.T 



(91,92,93)- 



The factorization of the divergent parts of the various virtual contributions, as multi- 
ples of the corresponding Feynman amplitudes at Born level, M$ , implies that the overall 
infrared and collinear divergences multiply the complete Born amplitude, M.b = Hi-M-^B ■ 
We can summarize our results for the virtual corrections to the individual fermion lines by 
writing the complete virtual amplitude M.y as 
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C F e 3 53 M$ V2 y 3>T g 2 , g 3 ) g? h g?*g?* + 0(e) 
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where the sums run over the different orderings of the attached weak bosons and the relevant 
topologies of Fig. ^ when using effective polarization vectors for the electroweak bosons, as 
discussed below Eq. (J2.5J) . Note that M.y is completely finite. The NLO contribution to the 
cross section at order a s comes from the interference of the virtual amplitude with the Born 
term. For corrections to a quark line it is given by 



2Re[M v M* B ] = \M B \ 



2vr 



-c. 



47T/1 



Q 2 



r(i + e) 



e z e 



+ 2 Re 



M V M* B 



(2.9) 

The divergent piece appears as a multiple of the Born amplitude squared and it cancels 
explicitly against the phase-space integral of the dipole terms (see Ref. [THj and Eq. (2.10) 
of Ref. 0) 



<I{e)>=\M B ?^^C ^ 
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which absorbs the real-emission singularities which are left after factorization of the parton 
distribution functions. After this cancellation, all the remaining integrals are finite and can, 
hence, be evaluated in d = 4 dimensions. 



2.4 Technical details 



Our Monte Carlo program computes all amplitudes numerically, using the helicity technique 
and the formalism of Ref. ^H]- For the tree-level and real-emission amplitudes (including 
counter-terms), the method is straightforward, since these contributions are finite at each 
phase-space point. The evaluation of the helicity amplitudes is very fast, due to the modular 
structure that one achieves by grouping the whole set of diagrams according to the topolo- 
gies illustrated in Fig. [T] The W + ^e + p e and W'^fi'v^ decay amplitudes and the single 
index leptonic tensors Ty (V = Z, 7) are effective polarization vectors which only depend 
on the lepton momenta. They are the same for all subprocesses, i.e. they do not depend 
on quark flavor or whether quarks and/or anti-quarks scatter. Similarly the second-rank 
leptonic tensors Ty V , T^ ±v and T^ +w _ are independent of quark flavor and come in just 
two kinematic configurations, depending on whether or not an external gluon is attached 
to the upper or the lower quark line in Fig. ^ Correspondingly, the leptonic tensors are 
calculated first in our numerical program and then used in crossed subprocesses and sub- 
traction terms. The code for these leptonic tensors has been generated with MadGraph and 
adapted to the tensor structure required for our full program. We note, in passing, that 
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this approach allows for straightforward inclusion of new physics effects in the electroweak 
sector: only the leptonic tensors would be affected by modifications like anomalous three- 
or four-gauge-boson couplings or strong electroweak-boson scattering. One major advantage 
of the modular strategy is the increase in computational speed. In the calculation of the 
real-emission contributions, which constitute the most CPU-time intensive part of the code, 
our program is about 70 times faster than a direct use of MadGraph-generated routines for 
the individual subprocesses. 

Special care has to be taken in the extraction of the finite parts M.ViV 2 ,t (qi, Q2) and 
M-ViV2V3,t (Qi, 12, 93), which are contained in the full virtual amplitude of Eq. ()2.8j) . In order 
to keep the expressions small and fast to evaluate, we have implemented the PV tensor 
reduction numerically. Here we are adopting a natural extension of the PV notation, and we 
call Eij the coefficient functions from the tensor reduction of pentagon integrals. Since the 
finite A4v!V2,r (<?i, 92) and M.ViV2V 3 ,t (<Zi, <?2, Q3) virtual sub-amplitudes only contain the finite 
pieces of the various tensor integrals, one needs to track how the divergent contributions in 
the expressions of the scalar integrals feed into the expressions of the tensor coefficients Bij, 
Cij, Dij and E^, and how they generate finite contributions in coefficients that contain a 
factor (d — 4) in the numerator. The resulting analytical expression for M. VlV2 T (^1,^2), in 
terms of finite functions, is given in Ref . |B] • We postpone to a future paper [22] any further 
technical discussion about the computation of M.ViV2V a ,t (Qi,Q2,Q3)- 

3 Checks and implementation in a parton-level Monte Carlo 

The cross-section contributions discussed in the previous section have been implemented 
in a fully-flexible parton-level Monte Carlo, which is very similar to the programs for Hjj 
and Vjj production in VBF as described in Refs. [7] and [S]. The matrix-element calcu- 
lation is divided into three main parts, that deal with the evaluation of the tree-level, the 
real-emission and the virtual contributions. All elements have been extensively tested as 
detailed below. 

Tree-level contribution 

We have compared our tree-level code with purely MadGraph generated output, and we have 
found agreement with a typical relative accuracy of 1CT 10 . 

Real-radiation contribution 

The same comparison has been performed for the real radiation contributions, with typical 
agreement at the 10~ 10 level. In addition, we have also checked the QCD gauge invariance of 
the real-emission corrections. More specifically, the real-emission amplitude for the process 
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qq'^qq'g e + z/ e /z has the form 

M R = e^p)M Rj (3.1) 

where p is the momentum of the emitted gluon and its polarization vector. Gauge in- 
variance demands that the amplitude A4r remains unchanged upon the substitution e^(p) — ► 
£fi(p) + Pp^ (with P arbitrary), that is 

Pfl M R = 0. (3.2) 
This relation is satisfied within the numerical accuracy of the program. 

Virtual contribution: code checks 

As far as the virtual contribution is concerned, we have implemented two different codes, 
one analytical, in MAPLE, and one numerical, in f ortran. The analytical code sums all the 
eight Feynman diagrams in Fig. which we call P M1At2Ai3 for uncontracted polarization vec- 
tors of the three electroweak bosons, and writes it in terms of the PV coefficient functions, 
Bij, . . . Eij, in d dimensions. We schematically represent this tensor reduction by 

V m (h, qi , q 2 , q 3 ) = £ r^ 2M3 , (3-3) 

u 

where (PV)^- = {B^, Cy, D^, E^} is one of the Passarino-Veltman coefficient functions, and 
the (finite) tensors correspond to spinor products describing the quark lines in Fig. El 

The Vmnzna and the (PV)ij still contain divergent contributions. We denote their finite parts 
by ^uiwzMsi Cij, Dij, E^, respectively. 

The analytic code contains all the recursion relations that can be used to reduce the 
PV coefficient functions to combinations of scalar integrals only: B , Cq, D and E func- 
tions. The Eq function can be further expressed in terms of the sum of five D Q functions, 
as described in Ref. j2H]; when d = 4 — 2e, in the limit e^O. The analytic continuation of 
Do functions was checked against Ref. [21]. The tensor reduction down to scalar integrals, 
and the direct substitution of the corresponding expressions computed in d = 4 — 2e dimen- 
sions have been used to check the structure of the divergent terms, and to show that, once 
contracted with the Born amplitude, they are given by Eq. (|2.8J) . 

The expression of V^^^ in terms of PV coefficient functions is turned by MAPLE into a 
f ortran code, where care is taken to obtain the correct limit when d— >A. All the technical 
details about this part of the program will be given in Ref. [22J. 

Both the analytical and the f ortran code have been checked extensively using gauge 
invariance, applied at different levels of complication: 

- at the level of the single pentagon loop (diagram (a) of Fig. EJ), 
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at the level of the sum of all the virtual corrections along a single quark line, "P m/l2At3 , 
and at the level of the entire scattering process (for the f ortran code). 



To illustrate an example of gauge check, we consider the simpler case of the pentagon loop 
of diagram (a) in Fig. |U1 



£MlM2M3(^l> 9i, 92, 93) 



d d l 



77 



1 1 

K23 /+ $1 " 



1 



(3.4) 

where g i2 = gi + 92, 9123 = 91+92 + 93, k 2 = ^123 + h- Gauge invariance is simply the 
statement that, upon contracting any one tensor index with the corresponding momentum, 
and expressing the contracted gamma matrix as the difference of the two adjacent fermionic 
propagators, the pentagon can be reduced to box integrals (see diagram (b) in Fig. [HJ) 

9i 1 £/^2M 3 ( /i + 9i,92,9s) = ^Ws ( k i, 9i + 92, 93) -£W 3 ( fc i + 9i,92,9s), (3.5) 
92 2 £ / ^2M3( A + 9i,92,9s) = ^Wa^i) 9i,92 + 93) -^1^(^1)91 + 92,93), (3.6) 
93 3 ^i^2M 3 ( fc i>9i,92,93) = V^^kx, q 1} q 2 ) -T>„ lfi2 (k 1 ,q 1 ,q 2 + q 3 ) ■ (3.7) 

Using the PV tensor reduction, we can express £ WM2At3 as a sum of coefficient functions up 
to (see Eq. fJ3.3|0 . and V^ u as sum of coefficient functions up to Dij, and generate the 
corresponding f ortran code for their finite parts. We can then check that the analytic or 
numeric expression for £ w/J2At3 , once an external index is contracted with the corresponding 
momentum, agrees with the right-hand-sides of Eqs. (|3.5J1 - (|3.7|) . Analogous relations hold at 
the level of the 7 ? Atl ^ 2//3 which represent the sum of all the virtual corrections along a single 
quark line. Both tests are a strong check on the correctness of the entire code. 

Finally, we have implemented two independent codes to compute the virtual corrections 
for the neutral-current contributions. The relative amplitudes agree within the numerical 
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precision of the two f ortran programs. 



Virtual contribution: numerical stability 

Gauge invariance has been used not only to check the entire code but is used every time 
that a virtual contribution is computed at a given point in phase space. When the dia- 
grams of Fig. EJare contracted with the leptonic currents which represent the W + ^e + u e and 
W^^fi^u^ decay amplitudes, the helicity amplitude has the generic form 

Jf 1 J^ 2 T '/^atts • (3-8) 

For example, in the computation of the virtual corrections for the sub-amplitude (a) in Fig.d 
we need to evaluate 

J+J-V^^ik^q+^-^o) , (3.9) 

where J + is the electronic current from the decay of a W + with incoming momentum q + , J_ 
is the muonic current from the decay of a W~ with incoming momentum q_ and qo is the 
incoming momentum of the neutral vector boson. We evaluate this expression by projecting 
the four-vectors J± on the respective momenta 

J£=x ± qt + r» ± , (3.10) 

in such a way that, in the center-of-mass system of the W + W~ pair, the vectors r± have 
zero time component 

r ± -(q + + q.) =0, (3.11) 

so that 

j±- («+ + »-) (312) 

Equation (|3.9jl then becomes 

J+ 1 J-V^^ih, q + , q-, q ) = 7+ 1 r^ 2 P WM2M (&i , q+, g_, q ) + box contributions, (3.13) 

where we have used Eqs. (|3.5|l - (|3.7jl . 

The projections of Eqs. (|3.1()j) - (|3.13jl reduce the magnitude of the coefficients multi- 
plying the pentagon loops and their overall contribution to the virtual corrections. This 
"true" pentagon contribution to the cross section, defined by the interference of the residual 
r w r A»2 "p^ i/J2 ^ 3 type terms with the Born amplitude, is called 05 below. Minimizing it is im- 
portant in view of the fact that, in the tensor-reduction procedure a la Passarino-Veltman, 
Gram determinants appear in the denominators of the PV coefficient functions. There are 
points in phase space where these determinants become small and numerical results become 
unstable. We have developed a strategy to interpolate over these critical points (see Ref. [221) 
and, to make sure that the numerical accuracy is not spoilt, we check, numerically, that the 
analogs of Eqs. (|3.5j) - (|3.7j) are satisfied for the tensors V^^^. In Fig. (a) we show the 
fraction of subprocess events, / (counting pentagon corrections to the upper and the lower 
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6 6 

Figure 7: Dependence of the pentagon contribution to the cross section on the 
maximal relative numerical error, 5, which is allowed for the Ward identities of the 
residual pentagon diagrams. The fraction f of subprocess events, where a pentagon 
subamplitude is discarded for numerical reasons, is shown in panel (a). Panel (b) 
gives the pentagon contribution, o~§ } to the cross section. Error bars correspond to 
the statistical error of the Monte Carlo integration. For comparison, the statistical 
error on the overall cross section for a high statistics run with ± 0.06% accuracy is 
indicated by the horizontal band. 



quark line as different subprocesses), where these Ward identities are violated by more than 
a fraction 5. Numerical instabilities of 5 = 10% or more, for example, affect about / = 13% 
of the generated subprocess events. For subprocess events with Ward identity violations ex- 
ceeding 5 we discard the numerically unreliable MviV 2 v 3 an d correct the remaining pentagon 
contributions by a global factor, 1/(1 — /), in order to compensate for the loss. As demon- 
strated in Fig. [7| (b), this procedure leads to a constant overall pentagon contribution, cr 5 , 
when varying S between 0.001 and 1000. Numerical instabilities become large for S > 10 4 . 
For our Monte Carlo runs we choose 5 = 0.1. Since the pentagon contribution amounts 
to less than 0.5% of the cross section for EW W + W~jj production with VBF cuts, the 
error that is induced by this approximation affects our final NLO results at an insignificant 
level. For comparison the shaded horizontal band shows the size of the overall cross-section 
uncertainty for a high-statistics run, with an overall statistical error of 0.06%. 

The tensor reduction of the box-type virtual contributions is quite stable, numerically. 
We have checked that the corresponding Ward identities, derived in a similar way as for the 
pentagons, are violated at only one out of 10 6 phase-space points by more than l%c. 
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The box- and pentagon-type virtual corrections, the finite M.v x v 2 an d M-v^Vz terms in 
Eqs. ()2.8|) and ()2.9|) . whose evaluation is cumbersome and time consuming, amount to less 
than one percent of the full cross section. Therefore, the statistical error of these contribu- 
tions affects the accuracy of the full result only marginally and the number of Monte Carlo 
events for the computation of the box and pentagon corrections can be reduced substan- 
tially with respect to the Born cross section and the leading c v i rt \M-B\ 2 virtual contribution 
in Eq. ()2.9|) : in our program, the Monte Carlo statistics is reduced by a factor 16 for the 
generic box contributions and by a factor 128 for the pentagon contributions after the projec- 
tions of Eqs. (|3.1U|) - (j3.13|) . These elements, together with the efficient handling of leptonic 
tensors and the other speed-up "tricks" described in the previous section, yield a fast code, 
which allows us to perform high-statistics runs with small relative errors on the full NLO 
cross sections and distributions. For example, it took about five days of CPU time on a 
3 GHz Pentium 4 PC to obtain an accuracy of l%o on the distributions shown in the next 
section. 

As discussed in detail in Ref. care has to be taken in the treatment of finite-width 
effects in massive vector-boson propagators. In order to handle diagrams where vector bosons 
decay, like W(pe + p v ) — > £(pe) + ve(p u ), a finite vector-boson width Ty has to be introduced 
in the resonant poles of each s-channel vector-boson propagator. However, in the presence 
of single- and non-resonant graphs, like (a) and (b) in Fig. [21 this introduces violations 
of electroweak gauge invariance in a sub-class of diagrams, which would hold in the zero- 
width approximation. In the past, different methods, such as the overall factor scheme [213 
and the complex-mass scheme [21)] > have been applied to overcome these problems. We 
resort to a modified version of the complex-mass scheme, which already has been used in 
Ref. [H]. We globally replace m v with m v — imyTy, while keeping a real value for sin 2 9w- 
This prescription has the advantage of preserving the electromagnetic Ward identity which 
relates the tree- level triple gauge-boson vertex and the inverse W propagator [2Zj. It thereby 
avoids large contributions from gauge-invariance-violating terms. 

Throughout the calculation, fermion masses are set to zero, because observation of either 
leptons or (light) quarks in a hadron-collider environment requires large transverse momenta 
and hence sizable scattering angles and relativistic energies. For consistency, external b- and 
t-quark contributions are excluded. 

We have used a diagonal form (equal to the identity matrix) for the Cabibbo-Kobayashi- 
Maskawa matrix, Vckm- This approximation is not a limitation of our calculation. As long 
as no final-state quark flavor is tagged (no c tagging is done, for example), the sum over all 
flavors, using the exact Vckm, is equivalent to our results, due to the unitarity of the Vckm 
matrix. 

The VBF cuts, discussed in Sec. HI force the LO differential cross section for W + W~jj 
to be finite, since they require two well-separated jets of finite transverse momentum. For 
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the NLO contributions, initial-state singularities, due to collinear q—^qg and g — > qq split- 
ting, are factorized into the respective quark and gluon distribution functions of the proton. 
An additional divergence is encountered in those real-emission diagrams, where a t-channel 
photon of low virtuality is exchanged, thereby giving rise to a collinear q^q^ singularity. 
We avoid it by imposing a cut on the virtuality of the photon, Q^ min = 4 GeV 2 . Events that 
do not pass this cut are considered to be part of the QCD corrections to the ^7 — > W + W~jj 
cross section, that we do not calculate here. 

For the computation of cross sections and distributions presented in the following section, 
we have adopted the CTEQ6M parton distributions with a s (mz) = 0.118 at NLO, and the 
CTEQ6L1 set for the LO calculation [2H]- The CTEQ6 parton distributions include b quarks 
as active flavors. However, since in our calculation all fermion masses are neglected, we have 
disregarded external b- and t-quark contributions throughout. 

We have chosen M z = 91.188 GeV, M w = 80.423 GeV and G F = 1.166 x lO^VGeV 2 
as electroweak input parameters. The other parameters, oqed = 1/132.54 and sin 2 ^ = 
0.22217, are computed thereof via LO electroweak relations. Final-state partons are recom- 
bined into jets according to the Ut algorithm [29J, as described in Ref. [SO], with resolution 
parameter D = 0.8. 

4 Results for the LHC 

The parton-level Monte Carlo program described in the previous section has been used 
to determine the size of the NLO QCD corrections to the EW W + W~jj cross sections at 
the LHC. Using the fcr algorithm, we calculate the partonic cross sections for events with at 
least two hard jets, which are required to have 

p Ti >20GeV, |%-|<4.5. (4.1) 

Here jjj denotes the rapidity of the (massive) jet momentum which is reconstructed as the 
four- vector sum of massless partons of pseudorapidity \rj\ < 5. The two reconstructed jets 
of highest transverse momentum are called "tagging jets". At LO, they are identified with 
the final-state quarks which are characteristic for vector-boson fusion processes. 

We consider the specific leptonic final state e + v e pTv^. One obtains the cross sections for 
the phenomenologically more interesting final state containing any combination of electrons 
or muons {e + e~vv 1 fi + fi~uj?, e^^vv, but neglecting identical lepton interference and ZZ 
final states) by multiplying our cross sections by a factor of 4. In order to ensure that the 
charged leptons are well observable, we impose the lepton cuts 

Pti > 20 GeV , \r] e \< 2.5 , AR j£ > 0.4 , (4.2) 
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where ARje denotes the jet-lepton separation in the rapidity-azimuthal angle plane. In 
addition, the charged leptons are required to fall between the rapidities of the two tagging 
jets 

yj,min ^ ^ Vj,max • v 4 "^) 

Backgrounds to VBF are significantly suppressed by requiring a large rapidity separation of 
the two tagging jets. We here impose the cut 

A Vjj = \Vh ~Vh\ > 4 • ( 4 - 4 ) 
Furthermore, we require the two tagging jets to reside in opposite detector hemispheres, 

Vh x Vh < > ( 4 - 5 ) 

with an invariant mass 

Mjj > 600 GeV . (4.6) 

The resulting total cross section receives two major contributions, arising from the Higgs 
resonance, via H^WW decays, and from the WW continuum, which effectively starts at 
the H^-pair threshold. Already for Higgs boson masses as low as 120 GeV, the resonance 
contribution is quite noticeable and, because of the strong dependence on m H of the H^WW 
branching ratio, this resonance contribution is strongly dependent on the Higgs mass. When 
trying to show results for the WW continuum only, we therefore impose the additional 
requirement 

m ww = yJipe+P^+Pu e +P^) 2 > m H + 10 GeV , (4.7) 

i.e. the four-lepton invariant mass must be above the Higgs resonance. The resulting cross 
section is representative of the continuum above any light Higgs boson resonance (m# below 
the IV-pair threshold). 

The scale dependence of the total continuum cross section, for a Higgs boson mass of 
rriH = 120 GeV, is shown in Fig. |HJ This figure shows the scale dependence of the LO and 
NLO cross sections, for renormalization and factorization scales, /ir and fip, which are tied 
to the W mass 

= fa m w , ^f = £f m w . (4.8) 

The LO cross section only depends on the factorization scale. At NLO we show three cases: 
(a) C,f = £r = £ (solid red line); (b) variation of the factorization scale only, £p = £, £r = 1 
(dot-dashed blue line); and (c) variation of the renormalization scale only £ R = £, £ F = 1 
(dashed green line). The NLO cross sections are quite insensitive to scale variations: allowing 
a factor 2 variation in either directions, i.e. considering the range 0.5 < £ < 2, the NLO 
cross section changes by less than 2% in all cases. Compared to this small variation, the 
factorization scale dependence of the LO cross section is quite sizable, amounting to a ±10% 
shift for 0.5 < £ < 2. Note that for /j,p — mw the LO cross section is only very slightly larger 
than the more stable NLO result, yielding a K factor K = uni^o/clo — 0.98, i.e. /ip = mw 
is an excellent choice for a LO estimate of the total continuum cross section. 
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Figure 8: Scale dependence of the total jj e + v e ^jL~i> fjL cross section at leading and 
next-to-leading order within the cuts of Eqs. \4- 1\ )- $71\ ) forpp collisions at the LHC. 
The contribution from the Higgs resonance (taken as mu = 120 GeV) is excluded. 
The factorization scale \xp and/ or the renormalization scale are taken as mul- 
tiples of the W mass, £mw, and £ is varied in the range 0.1 < £ < 10. The NLO 
curves are for fip = fJ-R = £,mw (solid red line), = f^w an d [ir = £ mw (dashed 
green line) and hr = mw and hf variable (dot- dashed blue line). The dotted black 
curve shows the dependence of the LO cross section on the factorization scale. At 
this order, there is no dependence on a s (fiR). 



Also for larger Higgs boson masses, m H > 2m w , the reduction of the scale dependence 
at NLO is comparable to the light Higgs case. However, since the resonance contribution 
can no longer be trivially separated from the WW continuum, we now show, in Fig. El the 
total cross section within the cuts of Eqs. (|4.1|) - (j4.6|) as a function of m# and for different 
scale choices, /i = ^rnw with £ = 0.5, 1 and 2. At NLO, the scale dependence is hardly 
visible while at LO one again finds a sizable factorization scale dependence. 

The small scale dependence which is observed for the total cross section at NLO is also 
found for infrared-safe distributions. Typically, scale variations between 0.5 mw and 2mw 
change distributions by about 2%, with somewhat larger variations, up to 6%, sometimes 
occurring in the tails of the distributions shown below. 

The K factor close to unity, which was found for the total cross section, no longer persists 
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Figure 9: Higgs mass dependence of the total pp^ jj e + v e cross section at LO 
and NLO within the cuts of Eqs. \4-l\) ~ \4-@il - Results are shown for renormalization 
and factorization scales [i = 0.5 mw , mw and 2m\y- 



for distributions. We demonstrate this effect by showing a few experimentally relevant 
distributions together with the dynamic K factor which is defined as 

= da NLO /dx 
do~LO I dx 

In the following the Higgs boson mass is taken as = 120 GeV and we show cross sections 
for the continuum above mww — 130 GeV and within the cuts of Eqs. (|4.1|) - (j4.7|) . All panels 
are for the scale choice fiF = = mw- 

A fairly strong shape change in going from LO to NLO is found for the tagging-jet 
transverse- momentum distributions. This is shown in Figs. El and ^2 where the larger and 
the smaller of the two tagging-jet transverse momenta are shown at LO (dashed black curves) 
and at NLO QCD (solid red lines), together with their ratio, the K factor of Eq. (J4.9j) . In 
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Figure 10: Transverse-momentum distribution of the highest-p? tagging jet in EW 
W + W~jj production at the LHC. In panel (a) the NLO result (solid red line) and 
the LO curve (dashed black line) are shown. Their ratio, the K factor as defined in 
Eq. 14- is shown in panel (b). 



particular the former, da / dp?^ , shows a clear shift to smaller transverse momenta at NLO, 
which corresponds to a K factor varying between 1.2 and 0.8 as p™tag increases from 20 GeV 
to 400 GeV. The effect for dcr/dp^ , in Fig. HH is slightly smaller, but still pronounced. 
The change in the jet transverse-momentum distribution also feeds into the shape of the 
lepton transverse-momentum distributions. In Fig. El we depict the transverse momentum 
for the hardest of the two charged leptons. Again small transverse momenta are enhanced 
at NLO, leading to a K factor between 1.04 and 0.84. 

In contrast to the transverse momentum distributions, angular distributions of the leptons 
are hardly affected by the NLO corrections. As an example, we show the azimuthal angle 
between the two charged leptons in Fig. The K factor is almost constant and equal to 
0.98. The typically large angle between the charged leptons is important for the reduction 
of W + W~jj continuum events in the search for H^-WW^l + l~$ T decays jlj. 

Another distribution which is important for the Higgs search at the LHC is the transverse 
mass of the l + l~vv system, which is defined as [3] 

Mjf w = fi$ T + E TM y - (p Ti , + j* T )2 , (4.10) 

where the transverse energies are given by 
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Figure 11: Same as Fig. \1(A but for the smaller of the two tagging-jet transverse 
momenta. 
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Figure 12: Transverse-momentum distributions of the hardest final-state lepton in 
EW W + W~ jj production at the LHC. In panel (a) the NLO result (solid red line) 
and the LO curve (dashed black line) are shown. Their ratio, the K factor as defined 
in Eq. is shown in panel (b). 
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Figure 13: Azimuthal angle separation between the two charged leptons for contin- 
uum jje + v e \i~Vn, events at the LHC. Curves are as in Fia. MIA 

#t = \l& + m lu « y/& + < ■ (4-11) 

While the invariant mass of the W + W~ pair cannot be reconstructed, due to the presence of 
neutrinos, Mjf w is fully accessible. The effect of NLO QCD corrections on Mjf w is again 
modest, as can be seen in Fig. HU 

5 Conclusions 

Vector-boson fusion at the LHC represents a class of electroweak processes which are 
under excellent control perturbatively. This has been known for some time for the most 
interesting process in this class: Higgs boson production via VBF has a modest K factor 
of about 1.05 for the inclusive production cross section |3T] and this result also holds when 
applying realistic acceptance cuts [Zj. Similar results were also found for Wjj and Zjj 
production via VBF [8]. 

In the present paper, we have extended these calculations to the electroweak process 
pp^e + h , e fi~u fJ- jj at NLO in QCD, when the final-state particles are in a kinematic config- 
uration typical of VBF events. This corresponds to leptonic final states in the vector-boson 
scattering processes VV~^W + W~ (V is a 7 or a Z) and W + W~—>W + W~, but with full 
NLO QCD simulation of the associated tagging jets. The calculation has been implemented 
in the form of a fully-flexible parton level Monte Carlo program and, thus, allows to imple- 
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Figure 14: Transverse mass distribution for the e + v e /j, Vn system in W + W jj 
events at the LHC. Curves are as in Fig. \M The definition of M^ w is given in 
Eq. 



merit completely general experimental cuts. The size of the QCD corrections is similar to 
those found for Hjj and Vjj production in VBF, and corresponds to a shift of a few percent 
in typical integrated cross sections expected for VBF cuts. Some distributions, however, are 
affected somewhat more strongly, with dynamical K factors ranging between 0.8 and 1.2, in 
particular for transverse-momentum distributions. At least as important is the stability of 
the NLO result: the residual scale dependence is at the 2% level for cross sections integrated 
within VBF cuts. 

The numerical code is quite fast, reaching permille level statistics on distributions within 
5 days of running on a standard 3 GHz PC. A l%o error on integrated cross sections is reached 
in about 1 day. This high speed has been obtained by avoiding the recalculation of recurring 
subamplitudes in different sub-processes contributing at a given phase-space point. A key 
ingredient is a modular structure of the numerical amplitude calculation which separates the 
weak-boson scattering sub-amplitudes into leptonic tensors, which can be changed without 
altering the validity of the QCD corrections. Such changes could reflect the inclusion of 
anomalous three- or four-vector-boson couplings or of any other new physics in weak-boson 
scattering. We leave such generalizations for the future. 
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